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We realized a quantum geometric “charge” pump for a Bose-Einstein condensate (BEC) in the 
lowest Bloch band of a novel bipartite magnetic lattice. Topological charge pumps in filled bands 
yield quantized pumping set by the global - topological - properties of the bands. In contrast, 
our geometric charge pump for a BEC occupying just a single crystal momentum state exibits 
non-quantized charge pumping set by local - geometrical - properties of the band structure. Like 
topological charge pumps, for each pump cycle we observed an overall displacement (here, not 
quantized) and a temporal modulation of the atomic wavepacket’s position in each unit cell, i.e., 
the polarization. 

PACS numbers: 


Ultracold atoms in optical lattices provide a unique set¬ 
ting for experimentally studying concepts that lie at the 
heart of theoretical condensed matter physics, but are 
out of reach of current condensed matter experiments. 
Here we focus on the connection between topology, ge¬ 
ometry, and adiabatic charge pumping GHZ] for Bose- 
Einstein condensates (BECs) in cyclically driven lattice 
potentials. 

Particles in periodic potentials form Bloch bands with 
energy e n (q) and eigenstates |T n (g)} = exp (iqx) \u n (q)) 
labeled by the crystal momentum q along with the band 
index n. The states \u n ) retain the underlying period¬ 
icity of the lattice, set by the unit cell size a. Motion 
in lattices is conventionally understood in terms of these 
bands: metals are materials with partially filled bands, 
while insulators have completely filled bands. In this 
context, a topological charge pump is a counterintuitive 
device, where charge motion - conduction - accompanies 
the adiabatic and cyclic drive of an insulating lattice’s pa¬ 
rameters. Thouless showed that this conduction is quan¬ 
tized, completely governed by the band-topology EE]. 
Although various charge pumps have been realized in 
condensed matter devices - such as modulated quan¬ 
tum dots [T0HT2] . ID channels driven by surface acoustic 
waves m, and superconducting qubits m - Thouless 
pumps remain unrealized in condensed matter settings 
but have been demonstrated in recent experiments with 
cold-atom insulators [T5J [16] . 

Here we break from this established paradigm for insu¬ 
lators and create a quantum charge pump for a BEC in 
a one dimensional (ID) lattice [TTUT91 occupying a sin¬ 
gle crystal momentum state q. This charge pump gives 
non-quantized motion sensitive to the Berry curvature at 
q integrated over the whole pump cycle, a local geomet¬ 
ric quantity, rather than a global topological quantity. 
Berry curvatures play an important role in condensed 
matter systems. An iconic example is the integer quan¬ 


tum Hall effect, where the electrons acquire an anoma¬ 
lous transverse velocity proportional to the Berry curva¬ 
ture and the quantized Hall conductance is given by the 
Berry curvature integrated over the whole 2D Brillouin 
zone (BZ) [20]; recent cold-atom experiments in 2D have 
measured such curvatures integrated over part [21.1 |22] 
or all 23 of the BZ. In an analogous way, ID lattice 
systems, driven cyclically in time £, have a generalized 
Berry curvature defined on the 2D effective BZ in q, t 
space. This curvature is the source of an anomalous ve¬ 
locity [24] , utilized to drive an adiabatic quantum pump¬ 
ing process. 

The Rice-Mele model [251428] of a bipartite lattice with 
a unit cell consisting of A and B sites is the paradigmatic 
system for understanding quantum pumps. The Hamil¬ 
tonian for this tight-binding model is 

Hrm = - ^2 [(t + St) bjCLj + (t - St) Oj +1 bj + H.c 

3 

+ AJ2 (a]aj - b]bj) , (1) 
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where d[- and describe the creation of a particle in 
unit cell j and sublattice site A or B respectively. The 
nominal tunneling strength t is staggered by St , and the 
sublattice sites are shifted in energy by A. 

We investigated quantum pumping in a novel ID 
(along e x ) bipartite magnetic lattice (building on 
Refs. [EH [30]) that in effect allowed independent control 
of £, St, and A. As shown in Fig. [lja)-(b), our magnetic 
lattice for 87 Rb arose from the interplay of one rf and 
two “Raman” fields that coupled the \f = 1; mp = ±1,0) 
“spin” states comprising the / = 1 ground state hyper- 
fine manifold, which were Zeeman split by hwz- The 
natural units of momentum and energy are given by the 
single photon recoil momentum hkR = 2i:h/\R and its 
corresponding energy Er = h 2 k 2 R /2m , where m is the 
atomic mass. In the frame rotating at the rf frequency 
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(a) Experimental geometry (b) Level diagram 




Position x/a rf phase (j )/tt 


FIG. 1: Bipartite magnetic lattice. (a,b) Dipole trapped 87 Rb 
BECs subject to a bias magnetic field Boe z had a Zeeman 
splitting o;z/27r — 0.817 MHz and a quadratic shift he = 
0.03 Er. These BECs were illuminated by four Raman beams 
and an rf magnetic field. Each of the two Raman couplings 
(strengths Q±) was derived from two cross-polarized Raman 
laser beams with frequency components uj and uj + Suj. (c) 
Adiabatic potentials colored according to (m x ) computed for 
h(n,n if ,5) = (6, 2.2,0 )E r , SQ/Cl = -0.1, and 0 = tt/ 4. The 
dashed curves plot the ±hCl x contributions to the potential 
experienced by states | m x = ±1). (d) Lowest two energy band 
energies plotted as a function of 0, otherwise with the same 
parameters as (c). 


Scj and under the rotating wave approximation, the com¬ 
bined rf/Raman coupling lead [31] to the overall Hamil¬ 
tonian 


H=^ + 0(x)-F + H q , ( 2 ) 

where F is the total angular momentum vector 
operator. We interpret Cl(x) = [Cl T f cos(</>) + 
Cl cos(2&r£), — fi r f sin(0) — SClsm(2kRx),y/2S]/\/2 as a 
spatially periodic effective Zeeman magnetic field, in 
which: fi r f is the rf coupling strength; Cl = + Cl_ and 

SCI m — Cl_ are derived from the individual Raman 
coupling strengths fi±; S = Suu — cuz is the detuning from 
Raman/rf resonance; <p is the relative phase between the 
rf and Raman fields. Additionally, Hq =—e(h 2 I—F 2 )/h 
describes the quadratic Zeeman shift, where I is the iden¬ 
tity operator. 

This spatially varying effective magnetic field pro¬ 
duces a ID bipartite lattice ElEl] with lattice constant 
a = Xr/2 with adiabatic (Born-Oppenheimer) potentials 
depicted in Fig.[l][c). This magnetic lattice is most easily 
conceptualized for small the Clcos(2kRx) term pro¬ 
vides periodic potentials for the \m x = ±1) states spa¬ 
tially displaced from each other by a/2 [dashed curves in 
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FIG. 2: Ground state spin projections, (a) Ground state spin 
projections at various <j) along with the predicted populations 
for h(Cl, SCl,Cl r f, S) = (4.4,0, 2.2, 0)Er. The associated adia¬ 
batic potentials [insets in (b)] have minima with spin projec¬ 
tion following the observed population’s trends, (b) Magne¬ 
tization derived from data in (a). 


Fig. B c)]; the resulting m x = ±1 sites are then staggered 
in energy, giving A « A max cos (0), with A max = fi rf / V2. 
The Cl y term couples these sublattices together: the 
rf term —Cl r { sin(0) generates constant height barriers 
(largely specifying £), which become staggered by the 
—SClsin(2kRx) contribution (largely specifying St). 

Figure [ljd) plots the energies of the resulting low¬ 
est two bands as a function of ej) (modulating A cosinu- 
soidally). Although our lattice is not in the tight binding 
limit, the band structure qualitatively matches that of 
the Rice-Mele model. In the remainder of this article, we 
focus on the lowest band n = 0 and will henceforth omit 
the band index. 

As illustrated by the shading in Fig.[ljc), in each unit 
cell the sublattice sites are “labeled” by their F x spin 
projection with \m x = — 1) site on the left and \m x = +1) 
site on the right. To confirm this, we adiabatically loaded 
| m z = 0) BECs into the lattice’s ground state by simul¬ 
taneously ramping the detuning from 5 Er to zero while 
ramping on the coupling fields in 10 ms. Following prepa¬ 
ration, our measurement sequence began with a n/2 spin 
rotation along e y , allowing us to measure the eigenstates 
of F x in our F z measurement basis. We achieved this tt/2 
rotation with a 44 fis pulse from an additional rf field 
with phase </> rot = n/2 and strength hCl r f >rot = 2.2E?#, 
applied while the Raman coupling was greatly reduced 
(Cl Cl r f ?ro t) and the lattice rf coupling was off (Cl r f = 0). 
We then abruptly removed the remaining control fields 
along with the confining potential and absorption imaged 
the resulting spin-resolved momentum distribution after 
a 20 ms time-of-flight (TOF) period in the presence of a 
magnetic field gradient along e y . 

Figure [ 2 ] shows the measured F x spin composition m 
and magnetization for adiabatically loaded BECs as a 
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function of 0 with SQ = 0. Because A(0) controls the 
relative depth of the \mn x = ± 1 ) wells, we observe ground 
states spin populations that follow this “tilt”. For exam¬ 
ple, when 0 = 0 or tt the double-well is strongly tilted 
and we observe the near perfect spin magnetization, con¬ 
sistent with atoms residing in the individual sub-lattices; 
in contrast, when 0 = 7 r/ 2 , the double-wells are balanced 
and we observe equal populations in each \m x ) state as 
expected for equal occupancy of both sub-lattices. Thus 
the magnetization [Fig. |2jb)] measures the mean atomic 
position within each unit cell, i.e., the polarization. 

Having constructed a physical realization of the Rice- 
Mele model, and demonstrated the requisite control and 
measurement tools, we now turn our attention to topolog¬ 
ical and geometrical charge pumping. These fundamen¬ 
tally quantum mechanical effects rely on the canonical 
commutation relation between position and momentum. 
Consider a finite wavepacket with center of mass position 
(COM) ( x) = (T|x|T), subject to a lattice Hamiltonian 
H that is adiabatically modulated with period T, i.e., 
H{t) = H(t-\-T). After one cycle, any initial crystal mo¬ 
mentum state is transformed |T(g)) —>> exp( 27 (g)) |T(g)), 
at most acquiring a phase, where q is the crystal momen¬ 
tum operator; this defines the single-period evolution op¬ 
erator U T = exp(ry(g)). The time-evolved position oper¬ 
ator UjixUt = x-dq^f(q) is displaced after a single pump 
cycle. 

The displacement is particularly simple in two limits: 
when just a single crystal momentum state is occupied or 
when every crystal momentum state in the BZ, —tt/ a < 
q < tt /a, is occupied with equal probability. As for our 
BEC, when a single |go) state is occupied the displace¬ 
ment is Ax(qo) = — d q ^(q)\ qo . Both the dynamical phase 
7d (g) = —c(q)T/h from the time-average energy e(g), 
and the geometric Berry phase 7 b( g) = i (u\dtu)dt 
contribute to 7 (g) = 70 (g) + 7 b(#)• In agreement with 
conventional descriptions HU M HU, this predicts a 
mean velocity v(q) = d q e(q)/h — T -1 F(q,t)dt. The 
first term is the usual group velocity and the second term 
- the anomalous velocity - derives from the Berry cur¬ 
vature F(g, t) = i({d q u\d t u) — ( d t u\d q u )). In our experi¬ 
ment, the BEC occupied the minimum of e(g, t) at g = 0 
during the whole pump cycle giving d q e(q) = 0 , so only 
the geometric phase 7b (g) contributed to the per-cycle 
displacement A x(q = 0) = — f Q F(q = 0, t)dt. 

In the contrasting case of a filled band, the aver¬ 
age group velocity is also zero and the displacement 
is Ax = —a f BZ 9g7B(g)dg/27r; this is often expressed 

as Ax = a J Q T 9t7z a k(^)dt/27r. The Zak phase 7 z a k = 
i f BZ (u\d q u)dq, a topological property of ID bands, is 
the Berry’s phase associated with traversing the ID BZ 
once, in the same way that 7 b (g) is a Berry’s phase taken 
over a pump cycle. 

Our lattice’s Zak phase is plotted in Fig. [3](a) ; this Zak 
phase is qualitatively indistinguishable from that of the 
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FIG. 3: Band geometry and topology computed for 
h(Q,Q T f,5) = (6, 2.2,0 )Er. (a,b) Zak phase and q — 0 Berry 
curvature showing the dependence on both 5Q/Cl and 0. In 
(b), the arrows show experimental charge pump trajectories 
in Fig. S b )- (c) Adiabatic potentials (displaced vertically 
for clarity) computed for a range of 0 constituting a com¬ 
plete pump cycle at SQ/Cl = —0.4 (left panel) and —0.8 (right 
panel). Filled circles mark the local energy minima. 


Rice-Mele model, with singularities at 0 = ±tt/2 and 
SQ = 0 , signaling topological phase transitions across 
these points. For filled band experiments, pumping tra¬ 
jectories encircling these points give quantized charge 
pumping fl5l[T6]. Figure ^b) shows the richly structured 
Berry curvature F(q = 0,0) relevant to our experiment, 
which will be explored next. 

For our charge pump experiments, we linearly ramped 
the pump control parameter 0(£) = 27rt/T, effectively 
modulating the lattice potential in two qualitatively dif¬ 
ferent regimes (separated by a critical |£D/D| ~ 0.63). 
In the first [Fig. |3][c), left panel] the sublattice sites rise 
and fall but the local potential minima are essentially 
fixed in space; in the second [Fig. [3](d) , right panel] each 
minimum is only present for part of the pump cycle (the 
potential appears to “slide” by Fa per cycle). As these 
schematics imply, the associated pumping process gives 
either no displacement, or a quantized per-cycle displace¬ 
ment =b a for classical trajectories [33] . In quantum sys¬ 
tems, however, geometrical pumping is controlled by the 
Berry curvature, giving non-quantized per-cycle displace¬ 
ments that can in principle take on any value. 
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(a) Magnetization evolution during charge pump 



(b) COM motion in charge pump (c) Displacement per cycle 



FIG. 4: Geometric charge pumping, (a) Magnetization measured while linearly ramping 0 with period T — 2 ms, along 
with the prediction for hCl m 6.38(2 )Er, MO = 4.50(2)Fr, and HCl r f = 2.20(3 )Er. (b) Displacement plotted versus 0/27T 
(number of pump cycles). Trajectories i to iii are taken at SCl/Cl = 0.7,0, and —0.7, respectively; in each case hCl ~ 6Er and 
hCl rf = 2.20(3 )Er. Solid curves: simulation of charge pump in the trap. The small displacement near 0 = 0 is introduced 
by our loading procedure, (c) Measured displacement Ax per pump cycle (symbols), along with the prediction obtained by 
integrating the Berry curvature over our pumping trajectory (solid curve). The uncertainty bars represent 95% confidence 
interval. 


We studied adiabatic charge pumping in this lattice in 
two ways: in the first we observed the F^-magnetization, 
giving the polarization within the unit cells, and in the 
second we directly measured the displacement Ax of our 
BEC. In both cases we loaded into the lattice’s ground 
state and linearly ramped 0 = 2tt t/T, driving the Hamil¬ 
tonian with period T 33]. As shown in Fig. Up), the 
magnetization oscillated with the T = 2 ms period, 
demonstrating the periodic modulation of polarization 
per cycle. In good agreement with our data, the solid 
curves in Fig. |4^a) show the predicted behavior given 
our known system parameters. This agreement persists 
to long times: for example after 50 pumping cycles (for 
t = 100 — 110 ms) the contrast is unchanged, confirming 
the adiabaticity of the process [33]. 

Lastly, we performed a charge pumping experiment 
by directly measuring the cloud’s position in-situ for a 
range of SO,/Cl. We obtained in-situ density distribu¬ 
tions using partial-transfer absorption imaging [34] in 
which ps 6.8 GHz microwave pulses transferred ~ 5% 
of the atoms from \f,m z ) = 11, — 1) to |2,0) where they 
were absorption imaged. This technique allowed us to re¬ 
peatedly measure the in-situ density distribution for each 
BEC. Each observed displacement was derived from dif¬ 
ferential measurements of the cloud position taken just 
before and just after the pumping process, rendering our 
observations insensitive to micron-level drift in the trap 
position between different realizations. 

Figure Qb) shows data taken for SCI/Cl = 0.7, 0, and 
—0.7 along trajectories i, ii, and iii, respectively, with 
both increasing and decreasing phases. Our data dis¬ 
plays two expected symmetry properties. First, since the 
displacement A x(q = 0) = — f F(q, 0)d0 depends on the 
sign of the acquired phase, the direction of motion is re¬ 
versed when the ramp direction is inverted. Secondly as 
shown in Fig. F(q = 0,0) is an odd function of 


SCI/Cl, so the direction of motion is also reversed when 
SCI/Cl —> —SCI/Cl. Thus Ax is an odd function of both 0 
and SCI/Cl, and as expected we observe no motion when 
SCI/Cl = 0. 

The displacement was markedly non-linear when the 
pumping time became comparable to our trap’s 80 ms 
period, showing the influence of the confining poten¬ 
tial [43^. We included the harmonic potential in our real- 
space simulations by directly solving the time-dependent 
Schrodinger equation for our system [35]. The simulated 
results [Fig. 0b), solid curves] agree with our observa¬ 
tions. To extract the per-cycle displacement due to geo¬ 
metric pumping, we fit the sinusoidal predictions of our 
model to each data trace, with only the overall ampli¬ 
tudes and a small vertical offset as free parameters, giving 
the short-time per-cycle displacement [33 . Figure [4][c) 
shows these per-cycle displacements for a range of Ra¬ 
man imbalances. 

The in-situ cloud typically had a Thomas-Fermi radius 
of 30 /im, corresponding to a small momentum width of 
0.004 kR for our BEC. We estimated the thermal fraction 
to be « 5% given by our « 20 nK temperature (momen¬ 
tum width of 0.24 &r). Moreover, the per-cycle displace¬ 
ment is nearly independent of q for \q\ < 0.25 kR [33]. 
These allow us to compare the data with the expected 
displacement from integrating q = 0 Berry curvature 
[Fig-0 c), solid line], showing an excellent agreement and 
confirming the geometric origin of our quantum charge 
pump. 

Our magnetic lattice enables new experiments with 
ID topological lattices. Berry curvatures at q ^ 0 can 
be probed by performing the charge pump pairwise at 
±\q\ (for example prepared via Bloch oscillations [36]). 
The dynamical phases in these cases are opposite and 
therefore cancel while Berry curvatures (even in q) con¬ 
tribute equally to the displacements [33]. Furthermore, 
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protected edges states, a hallmark of topological systems, 
will be present at the interface between regions charac¬ 
terized by different topological invariants j37H39j . Since 
in our lattice the topological index is set by the rf phase, 
a bulk topological junction can be generated by replacing 
the rf field with an additional co-propagating pair of Ra¬ 
man laser beams in which just one beam has an abrupt it 
phase shift in its center. This provides a static model of 
the soliton excitation mode in polyacetylene [25| l4Q]. Ter¬ 
minating our lattice with hard-wall boundaries will give 
rise to similar end states - somewhat analogous to Majo- 
rana fermions in ID topological superconductors [381 SB 
- with a spin character. 
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SUPPLEMENTARY MATERIAL 
RELEVANT RAMAN COUPLING 

The Raman coupling resulting from the vector light 
shift (proportional to E x E) is non-zero for our cross- 
polarized counter-propagating lasers, but is zero for the 
otherwise resonant contributions with each beam taken 
separately. For example, two Raman beams propagat¬ 
ing along +e x with frequencies uj and uj + Suj have the 
same polarization along e z in Fig.|5ja) and hence give no 
vector light shift. The same frequency components from 
the cross-polarized Raman beams [e.g., beams with fre¬ 
quency uo propagating along -\-e x with polarization along 
e z and propagating along —with polarization along 
e y in Fig. [5ja)] do not provide a time dependent vec¬ 
tor light shift for coupling the magnetic levels (or far 
off resonant from the Zeeman splitting in the rotating 
wave approximation). Therefore, the only two relevant 
Raman coupling pairs are marked in Fig.[5|a), with cou¬ 
pling strengths denoted as Q±. 

Figure [5|b) shows the level diagram of the three mp 
states in our 87 Rb BECs / = 1 ground state manifold as 
coupled by the experimental combination of Raman and 
rf coupling. For example, atoms starting at rest (ki = 0) 
in the |nrip = —1) spin state are coupled to the \mp = 0) 
state in three ways: the rf field provides coupling but 
leaves the momentum unchanged at hk = 0; while the 
two photon Raman transitions with strength Q± change 
the momentum to hk = ±2 hkR. 

We independently controlled the amplitudes of Cl± us¬ 
ing two rf sources with frequencies uj and uj + Suj, pro¬ 
duced from a Direct Digital Synthesized (DDS) signal 
generator (Novatech, Model 409B). We generated the de¬ 
sired Raman frequency pairs by splitting the output of 
each of the rf sources into two and controlled these four 
amplitudes individually. We identify these four rf signals 
by A+(uj), A-(uj), A+(uj + Suj), and A_(cj + Suj). Here, 
A+ and A_ denote the rf signal amplitudes. We then 
combined A+(uj) with A _(uj + Suj) and used the result¬ 
ing signal to drive the acousto-optic modulator for the 
Raman beams propagating along +e x with polarization 
along e z in Fig.[5|a). Similarly we combined A+(uj-\-Suj) 
and A_ (uj) to drive the acousto-optic modulator for the 
Raman beams traveling along — e^. 

LATTICE CHARACTERIZATION VIA PULSING 
EXPERIMENT 

We first characterized the magnetic lattice by study¬ 
ing the time-evolution of an initial \m z = 0) BEC follow¬ 
ing the abrupt turn-on of different combinations of our 
coupling fields, as shown in Fig. [6] After holding these 
fields constant for a time St we abruptly removed them 
along with the confining potential. We absorption im¬ 


aged the resulting spin-resolved momentum distribution 
after a 20 ms time-of-flight (TOF) period in the pres¬ 
ence of a magnetic field gradient along e y . Absent rf 
coupling we observed sinusoidal Rabi-like oscillations be¬ 
tween the \m z = 0) and \nn z = ±1) states (Fig. [6j left 
panel). While completing the lattice with the rf coupling 
makes the time-evolution more complex, it is still in good 
agreement with our predictions (Fig. [6| right panel). The 
individual coupling strengths were separately calibrated 
using a similar technique, and by fitting data in the right 
panel of Fig. [6] with as the only free parameter, we 
readily calibrated the rf phase. 

LIGHT-ATOM INTERACTION AND CONTROL 

OF RF PHASE 

The light-atom interaction due to a two-photon Raman 
transition can be expressed as 

Hi~a = &r,± • F. (3) 

The effective magnetic field from a single Raman coupling 
is 

SIr,± = 4=[fi± cos(2 k R x), -f2± sin(±2k R x), V25\, (4) 
v 2 

where the two coupling fields have an opposite sign on 
the e y component due to the opposite momentum kick. 
The rf magnetic field contributes the additional coupling 
f2 r f = l/\/2[Drf cos(0), — U r f sin(0), y/25\, where the rf 
phase (f) is defined relative to that of the Raman cou¬ 
pling. The effective magnetic field in Eq. (2) of the main 
text has the contributions from the above three coupling 
fields, and f. The DDS signal generator 

described in the previous section has four rf phase syn¬ 
chronous outputs, and we used two of the rf sources to 
generate frequency shifts on the Raman lasers. In ad¬ 
dition, we used the other two rf sources to generate the 
rf couplings for lattice and tt/2 spin rotation. The path 
lengths were different for the signal cables which drove 
the Raman AOMs and the rf coil. In addition, Raman 
beams after AOMs were sent to the chamber through 
optical fibers. The phase difference between the Raman 
couplings and rf coupling controlled by the DDS signal 
generator was different from the relative phase experi¬ 
enced by atoms. We therefore used both the pulsing and 
adiabatic loading measurements (Fig. [6] and Fig. 2 in the 
main text) to calibrate the relative phase on atoms at 
a certain commanded phase difference. After calibrat¬ 
ing the phase shift introduced by the different paths, we 
can add the known phase shift value to the DDS sig¬ 
nal generator. In our charge pumping experiment, we 
periodically modulated the lattice potential by linearly 
ramping the rf phase. To achieve this, the rf frequency 
was set to Suj r f = Suj A Suj' , corresponding to a linearly 
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(a) Experimental geometry 
e* 



(b) Level diagram 



FIG. 5: (a) Experimental geometry, (b) Level diagram of magnetic levels dressed by two Raman couplings with strength Q± 
and one rf coupling with strength Q r f. 
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FIG. 6: Time evolution of population in | m z — 0, ±1) states following the instantaneous turn-on of the lattice. Left panel: 
Simple evolution when only the Raman fields are applied with hQ+ = 5.45(2)F± and Ml- = 0.94(1 )Er. Right panel: Evolution 
when an rf field with strength Ml r f s= 2.20(3 )Er and phase 0 as — 0.02(1 )7r is added. In each case the solid curves depict the 
time evolution expected from our Hamiltonian using calibrated couplings; only 0 is a fit parameter in the right panel. 


varying phase relative to that of the Raman fields with 
T = 2 tt/Suj' . We switched between linearly increasing or 
decreasing phases by inverting the sign of Suj' . 


CLASSICAL, GEOMETRIC, AND 
TOPOLOGICAL PUMPS 

Our BECs had a typical temperature of 20 nK with 
~ 5% of the atoms in the thermal component. To show 
the contrast between a geometric pump and a classical 
pump, we first simulated trajectories of particles at an 
initial temperature of 0 K, which is the closest condition 
for describing the majority of atoms occupying the low¬ 
est energy state of the lattice. The classical trajectory is 
governed by the force derived from the lattice potential. 
Figure [7][a) shows the simulated per-cycle displacement 
versus for a geometric pump and a classical pump. 

We also included the topological pump result, which was 
obtained by integrating the Zak phase [Fig. 3(a) in the 


main text] over rf phase, describing the displacement 
of a filled band. The topological pump displays quan¬ 
tized per-cycle displacement as expected; the sign of the 
displacement depends on the sign of Raman imbalance, 
since the Zak phase in Fig. 3(a) of the main text also has 
mirror symmetry across zero Raman imbalance. 

The classical pump also displays quantized displace¬ 
ment. The transition point from zero to non-zero per- 
cycle displacement occurs at £U/U « ±0.63, which is 
consistent with the underlying motion of the adiabatic 
potential [Fig. 3(c) in the main text]. For large Raman 
imbalances, the adiabatic potential is displaced by one 
site per cycle, shifting the classical particle by one site as 
well; for small imbalances, two sublattices always main¬ 
tain local energy minima, giving no net force on the clas¬ 
sical particle released from local energy minima. 

We also computed the classical pump for atoms at non¬ 
zero temperatures: 30 nK and 100 nK. The initial ther¬ 
mal distribution was taken to occupy a single sublattice 
centered at the origin (which has the lower energy at 










(a) (b) 



5D/U SQ/Q 

FIG. 7: (a) Computed per-cycle displacement versus Raman imbalance SQ/Q for geometric, topological, and classical pumps. 
Classical particles have zero initial kinetic energy. Grey dashed lines indicate the critical Raman imbalance values SQ/Q = 
±0.63. (b) Computed per-cycle displacement versus 5U/U for classical ensembles at finite temperatures (30 nK and 100 nK) 
is shown in symbols. Ramping period for classical pumps is 2 ms. 


0 = 0 compared to the other sublattice regardless of the 
Raman imbalance), and the resulting per-cycle displace¬ 
ments are shown in Fig.^b). The 30 nK system is nearly 
quantized over a wide range of Raman imbalances, while 
the 100 nK ensemble deviates by more from the quantized 
values. In either case, the qualitative dependence on the 
Raman imbalance is in stark contrast to the geometric 
pump (solid line). 


BERRY CURVATURE AT NON-ZERO CRYSTAL 
MOMENTUM 


In our charge pump experiment, we prepared the BEC 
at q = 0 with a small momentum width of 0.004 k^. 
Therefore, we only consider the Berry curvature F(q = 
0,0) in the main text. BECs in a lattice at different 
crystal momenta can be prepared experimentally using 
Bloch oscillation |36| (such as applying optical potential 
gradient), prompting us to investigate the Berry curva¬ 
ture of our lattice at q ^ 0. Figure [ 8 ] plots the Berry 
curvature F(g, 0 ) at several q values versus the two ex¬ 
perimental parameters, Raman imbalance 5U/U and rf 
phase 0. To evaluate the Berry curvature, we used the 
discrete approach [18]. First, we defined the phase angle 
7 between two states as exp(—iq) = (u(R)\u(R + SR)) = 
1 ± 8R(u\8ru) : where R = q or 0. By taking the log¬ 
arithm of the previous equation, we obtained 7 = — Im 
\og(u(R)\u(R + 8R)) = 5Ri(u\dRu). Therefore, we can 
evaluate F(g, 0 ) as follows: 


F(q,<t>) = iddquldfu) - (d<f,u\d q u)) 

= i—{u\d<j>u) - i—(u\d q u) 

= ^^[-Imlog(u 2 |u 4 ) +Imlog(ui|u 3 )] 

- ^^[-Imlog(u 3 |u 4 ) +Imlog(ui|u 2 )], (5) 

where we define \ui) = \u(q, 0)), |iz 2 ) = 

\u(q + 8q, 0)) , \us) = |?i(g,0±50)), and \u 4 ) = 

\u(q ± 5q, 0 ± 50)); we used Sq = 10 _3 &r and 50 = 10 -3 
for computation. When q approaches the end of 
BZ, the largest Berry curvature value increases and 
F(q, 0) becomes highly concentrated near points with 
zero Raman imbalance and 0 = (2 n ± 1) x tt/2. 

Band gaps close at these points in the end of BZ, 
as shown in Fig. HRa). Topological phase transition 
occurs when these points are traversed where Zak 
phases possess singularities [Fig. 3(a) in the main 
text]. In addition, Berry curvatures are even in q , 
i.e. F(q, 0) = F(—q, 0). We want to show next that 
F(g,0) even in q is a general consequence of a ID 
lattice possessing a time-reversal symmetry. Note that 
time-reversal symmetry operates on the ” frozen” lattice 
at any moment during pumping (lattice Hamiltonian 
during a pump cycle does not satisfies H(—t ) = H(—t), 
which would give zero displacement). The Bloch 
states at q and — q can at most differ by a phase, i.e. 
Tu q ^{x) = u-q t <f,(x)* = exp(i/3(g, (t>))u q ^ (x), where f is 
the time reversal operator. Using this gauge and Eq. [5j 
we can show F{—q,(j)) = F(^,0) + \/{8<j>8q)[— (0 4 — 
02 ) + (03 - 0 i ) + (04 - 03 ) - (02 - 0 i )] = F(q, 0 ). Here, 
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(a) F(q = 0.25k R , (f>) 


(b) F(q = -0.25k R ,(f>) 



(c) F(q = 0.5 k R ,4>) 
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FIG. 8: Berry curvature F(q,<p) versus Raman imbalance SQ/Q and rf phase (p/n. (a-b): q = ±0.25/c#; (c-d): q = ±0.5&r; 
(e-f): q = ±0.75 k R - (g-h): q = ±0.9 k R . h(Cl, Or, S) = (6, 2.2, 0 )E R . 


exp (i/3j) is gauge term for state | Uj). This yields Berry 
curvatures which are even in q as well. 

The per-cycle displacement at different q values can 
be obtained by integrating the Berry curvature over the 
rf phase in a pump cycle A x(q) = — / F(q, (p)d(j), result¬ 
ing in Fig. [9j Since the Berry curvature is even in q, we 
only showed the displacement for positive q. The per- 
cycle displacements versus Raman imbalance are similar 
between q = 0 and q = 0.25 Icr: when reducing the im¬ 
balance value from = 1 —>• 0, the displacement 

value gradually reduces from less than \a\ to zero. At 
q = 0.5 kR, the displacements reach ±a at large Raman 
imbalances. For q > 0.5 kR, the displacement demon¬ 
strates a significantly different dependence on Raman 
imbalance: displacement value increases first when re¬ 


ducing the imbalance value. The displacement value is 
much larger than \a\ when \q\ approaches kR. This re¬ 
flects the large and highly concentrated Berry curvatures 
near zero Raman imbalance, shown in Fig. [8jg-f). 

When averaging the displacements of the geometric 
pump results for the entire BZ, we should obtain the 
result of a topological pump, as shown as black crosses 
in Fig. [7^a). Since the topological pump yields quan¬ 
tized displacement ±a depending only on the sign of Ra¬ 
man imbalance, the non-classically large displacement at 
low Raman imbalance and large q in Fig. 9] compensates 
the small displacement at low q. In Fig. lOj we showed 
the Berry curvature F(q, (p) on the plane spanned by 
crystal momentum q and rf phase 0, which effectively 
forms a torus for a filled band during charge pumping. 
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FIG. 9: Per-cycle displacement versus Raman imbalance 
for different crystal momentum values. Per-cycle displace¬ 
ments were obtained by integrating Berry curvatures in 
Fig. |8ja,c,e,g) over the rf phase. 


Berry curvatures are concentrated toward q = =b Jzr at 
< j) = ±7r/2, which is more obvious for low Raman im¬ 
balance [Fig. [To|[a-b)]. In addition, low Raman imbal¬ 
ance has a larger \F(q, (j))\max than high Raman imbal¬ 
ance. We integrated F(q, <fi) in Fig.[l0]to obtain the per- 
cycle displacement and indeed obtained quantized dis¬ 
placement of ±a(—a) for SCI /Cl < 0 (SCI /Cl > 0), which is 
consistent with the topological pump result [black crosses 
in Fig. 0a)]. 


BAND STRUCTURES 

The band structures of the lowest two bands are shown 
in Fig. [Tljfor three Raman imbalances. For each imbal¬ 
ance, we showed the band structures as a function of 
rf phase in a pump cycle (from top to down). We first 
noticed that the band structure is even in crystal mo¬ 
mentum at any </>, i.e. e(q) = e(—q). In general, a ID 
lattice Hamiltonian preserving time-reversal symmetry, 
such as the SSH model, yields a band structure which is 
symmetric in q. 

From the above two properties, our pumping scheme 
can be extended to probe Berry curvatures at q ^ 0 
from a pairwise measurement. Using Bloch oscillations 
[36], we can prepare BECs at any ±|go| and ramp the 
rf phase as we showed in the main text. Since the dy¬ 
namic phases from ±|g 0 | are opposite (symmetric band 
structure in q gives opposite group velocities) but the ge¬ 
ometric phase contributions are the same, the sum of the 
per-cycle displacements at ±|go| on ly probes the Berry 
curvature at qo. Specifically, Ax(+\qo\) ± Ax(—\qo\) = 
- f F(+\q o \,4>)d0- f F(-\q o \,4>)d0 = -2 f F(\q o \^)d0. 


After repeating the measurement for the several q values, 
one can obtain per-cycle displacements shown in Fig. [9| 
and hence map out the geometry of the entire band. 


FITTING PROCEDURE TO OBTAIN 
PER-CYCLE DISPLACEMENT AND STATE 
PREPARATION 


We expected the displacement to have a linear depen¬ 
dence on pump cycles from the geometric pumping with 
a per-cycle displacement given by A x(q = 0) = v x T = 
— f F(q , (j))d(j). However, the restoring force from the 
dipole trap would impact the actual displacements. We 
performed real-space simulations of the charge pumping 
process using the split operator method in the presence 
of a harmonic trapping potential, as shown in Fig. [12] 
for Raman imbalance ±0.7 (thick lines). To extract the 
contribution of geometric pumping, we can describe the 
simulated displacements as 

Ax± = v±/Vk± sin (Vk± x </>) + c±, (6) 


where subscripts ± represent the directions of increasing 
and decreasing phases. Here, k± are the effective spring 
constants from the trapping potential and magnetic lat¬ 
tice potential. 

During charge pumping, the phase was linearly ramped 
< j) = 2tt t/T. To linearly ramp the rf phase, we set a small 
frequency difference between the rf and Raman coupling 
fields as described before. Due to technical constraint 
on the DDS signal generator, we had to introduce the 
frequency difference in the beginning of the state prepa¬ 
ration. During state preparation, we ramped the cou¬ 
pling fields up to constant values and an initial detuning 
(Si = 2.5 Er) down to zero. We found through simula¬ 
tions that ramping the rf phase during state preparation 
introduced small offsets c± (Note that rf phases remained 
at constant values and Si = 5Fr during state preparation 
for Fig. 2 of the main text). Thin lines in Fig. 12 show 
the fitting results using Eq. [6] The difference between 
the average value of v± from the fits and the expected 
per-cycle displacement from integrating Berry curvature 
is < 1%. 

We can further fit our data at non-zero Raman imbal¬ 
ances [such as symbols in Fig. 4(b) of the main text] to 
the functional form: Ax = F +je xp/v / ^+ sin(\/fc + x 0) + 
c + ± C for 0 > 0 and Ax = F_ )exp /\/&- sin (Vk- X <t>) + 
C- ± C for < 0, where the three fitting parameters 
are F± ?exp and an overall offset C; the rest of the pa¬ 
rameters are taken from the fit results to the simulated 
displacements. The mean of F± )e x P extracts the per-cycle 
displacement from our data. For zero Raman imbalance, 
we simply fitted a linear line to the data shown as black 
symbols in Fig. 4(b) of the main text. 
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FIG. 10: Dependance of Berry curvature on q and 0 for four Raman imbalances. h(Q,Q T f,S) = (6, 2.2,0 )Er. 


FIELD STABILIZATION point. 


We took two in-situ partial-transfer images using two 
microwave pulses which were ±500 Hz detuned from ujz 
in all of our measurements. Based on the asymmetry 
(ratio of difference to sum) of two atomic signals, we ac¬ 
tively stabilized the laboratory bias magnetic field after 
each measurement. This can only correct a long term 
field drift in the laboratory. We rejected the data point 
if the asymmetry was larger than 0.25. Based on our 
independent field calibrations from pulsing a single Ra¬ 
man or rf coupling field alone, the post-selected atoms 
experienced a detuning which was within the quadratic 
Zeeman shift of he ~ h x 100 Hz. We intentionally in¬ 
troduced some detuning Si (by changing the bias field by 
a known amount) and then ramped it down during the 
state preparation, which can adiabatically load the BECs 
into the ground state of the magnetic lattice. 

For Raman imbalance of ±0.7 [trajectory i and iii in 
Fig. 4(b) of the main text], we took 30 images for each 
data point; for zero imbalance [trajectory ii in Fig. 4(b)] 
and imbalance of ±0.35 (trajectory not shown) we took 
10 images for each data point. Since we rejected data 
points based on the partial-transfer atom signals, we 
ended up keeping 60% — 70% of the images for each data 


ADIABATICITY OF CHARGE PUMP 


As shown in Fig. 11, the smallest gap at q = 0 is ~ 1 Er 
occurring at zero imbalance and (/> = (2n± 1) x 7 t/2. Our 
pumping speed of 500 Hz for non-zero Raman imbalance 
is appreciably slower than Er/H = 3.65 kHz, ensuring 
the adiabaticity. In addition, we had tried a factor of 
two larger ramping speed than 500 Hz in the experiment 
at Raman imbalance of ±0.7 and observed no significant 
difference in both the displacement and magnetization 
measurements. 


After each measurement of displacement, we also 
recorded the time of flight image of the same BEC. The 
contrast of the diffraction pattern (showing the distribu¬ 
tion in both momentum and spin) remained sharp dur¬ 
ing our pumping experiment: a good added indication of 
an adiabatic experiment. In addition, the magnetization 
measurement shown in the manuscript displays a good 
contrast for more than 50 pump cycles at the ramping 
rate we chose, confirming again the adiabaticity of our 
pumping process. 
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(a) 8 n/n = o (b) 5 n/n = -0.4 
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FIG. 11: Band structures of the lowest two bands of the magnetic lattice at three Raman imbalances: (a) SQ/Q = 0, (b) 
SQ/Q = —0.4, and (c) 5Q/Q = —0.8. For each Raman imbalance, we plot a series of rf phases <f> (from top to down) consisting 
of a pump cycle. Q r t , S) = (6, 2.2,0 )Er. 
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FIG. 12: Simulated displacements (thick lines) in charge 
pumping for Raman imbalance 5Q/Cl = ±0.7. h(Cl,Q T f,5) = 
(6.38, 2.2, 0 )Er. Thin lines show the fitting results of the sim¬ 
ulated displacement using Eq. [6] 



